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Recently a non-perturbative formula for the RG flow between UV and IR fixed 
points of the coefficient in the trace of the energy momentum tensor of the Euler 
density has been obtained for N = 1 SUSY gauge theories by relating the trace 
and R-current anomalies. This result is compared here with an earlier perturba- 
tion theory analysis based on a naturally defined metric on the space of couplings 
for general renormalisable quantum field theories. This approach is specialised to 
N = 1 supersymmetric theories and extended, using consistency arguments, to ob- 
tain the Euler coefficient at fixed points to 4-loops. The result agrees completely, 
to this order, with the exact formula. 
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The work of Seiberg [1] and many other authors has greatly clarified the non- 
perturbative structure of N = 1 supersymmetric theories. There is considerable evidence 
that there exist infrared attractive fixed points in many models, where beta functions van- 
ish and which define non-trivial superconformal invariant theories. As a consequence there 
has been renewed interest in the quest for a four-dimensional version of the Zamolodchikov 
c-theorem, which defines a function of the couplings for two-dimensional field theories which 
decreases monotonically under RG flow to the infra red limit and coincides at fixed points 
with the Virasoro central charge of the associated conformal theory. 

In a four-dimensional theory coupled to a background metric, the external trace anom- 
aly of the energy-momentum tensor is 

16tt 2 T a a = operator terms + c C apl6 C aPlS - a R a ^ s R* a ^ 5 , (1) 

where C a0l5 is the Weyl tensor and R* af3 lS = \e a p er) e^ S(T PR er ' < ap , so that R a ^ 5 R* a ^ s is the 
Euler density. Cardy has suggested that the Euler anomaly coefficient a(g) as a function 
of the running couplings g 1 (/i) provides the desired c- function. The total flow auv ~ a m 
between UV and IR fixed points has been calculated in both non-supersymmetric |J and 
more recently in many supersymmetric models ||,|5]|| and has been found to be positive in 
all models for flows from the trivial asymptotic freedom UV fixed point, while cuv ~ °ir 
has no definite sign in the models studied. 

In the analysis of [f||| the values auv and cuv were obtained from the free field 
content of the asymptotically free theories studied, while am and cjr were related to the 
R-current anomalies usually computed in studies of N = 1 duality. At a fixed point, where 
the renormalisation group beta functions vanish, by virtue of superconformal invariance 
the R-charges of chiral fields are related to their scale dimensions so that the result for 
the flow of a can be simply written in terms of the anomalous dimensions 7 r , at the 
IR fixed point, of the chiral superfields in the theory. With r denoting the gauge group 
representation R r , of dimension dimi? r , the non-perturbative formula, for all models with 
a unique R-current, is 

8(a uv - a IR ) = ^ dinLR r 7 r 2 (l - §7 r ) , (2) 

r 

(Analogous formulae for the flows of c and of flavour current anomalies are given in || 
but they are not relevant here). The formula (0) is a non-perturbative result, but it is a 
desirable check of the methodology to compare it with perturbative calculations which is 
the purpose of the present note. For models with a perturbatively accessible fixed point, 
the 7 r and also the coefficient a in ([D can be expressed as a power series in the couplings 
so that comparison can be made. 
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We now discuss the perturbative formalism of 0, in which the usual RG equations 
describing the variation of an arbitrary renormalisation scale (jl were extended to local 
rescalings of the background space-time metric. In this approach a can be calculated to 
high orders of perturbation theory in a general renormalisable quantum field theory, thus 
generalising earlier results in special cases [||. (The present notation for the anomalies 
in ([!]) is related to that of M by a = 167r 2 /3b, c = — 167r 2 /3 a .) The initial study [[7] used 
dimensional regularisation, but the results were also shown [Q to follow from consistency 
conditions related to the commutativity of local rescalings. In a renormalisable theory with 
couplings g 1 and associated beta functions (3 1 , it was shown how to define, to all orders in 
perturbation theory, a metric gjj = gjj on the space of couplings and associated one-form 
Wi so that 

d I 8a = t IJ p J , t I j = gu + d I Wj-djW I , 8a = 8a + Wi(3 I . (3) 

gu and also Wi may be calculated in terms of connected one particle irreducible vacuum 
graphs in the background geometry and (|||) then permits calculations of the perturbative 
expansion for a. The lowest order corrections to the one loop free field results occur in 
general at three loops and were given in [|7| . If Wj is an exact one- form then (^|) represents 
a gradient flow . Clearly at a fixed point d is equal to a and under renormalisation flow 
(H) implies 

±gi = -pigt) =► ^Sofo) = -guW 1 ^^) , (4) 

Thus, if the metric is positive, d(gt) monotonically decreases, as does the Zamolodchikov 
c-function in two dimensions. The positivity properties of the metric are not known in 
general but the leading perturbative contribution is positive, so one can establish the 
c-theorem at least in some region of weak coupling. It should be noted that it is the 
associated quantity d which is monotonic and not the anomaly a itself (which actually 
increases as one flows away from the UV fixed point in some models) . Further the formula 
(0) was derived under the condition that /3-functions vanish, so our comparison can only 
be made at an IR fixed point. 

In this note we apply the method of J7|,|8[ to general N = 1 supersymmetric theo- 
ries, assuming a simple gauge group G with coupling g and chiral G-invariant couplings 
Yijk-, Y % i k so that the superpotential is ^Y^kcj) 1 ^ (f) k for <\> % complex chiral superfields act- 
ing on which the generators of G are t a = —tj. In an appendix the 3- loop results of 
are briefly reviewed and specialised to the case where gauge, Yukawa, and quartic scalar 
couplings are related by supersymmetry. This gives a straightforward check of the formula 
(2) to 3-loop order. However the integrability conditions following from (|^) allow, with- 
out any new Feynman graph calculations, the expression for the metric to be extended to 
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higher order and then 4-loop results for a to be found. To achieve this we use the NSVZ 



beta function |TjJ for the gauge coupling g to three loop order. 

It is convenient to eliminate factors of 4n by writing g = g/Air, Yijk = Yijk/Air, Y^ = 
Y 1 ^ /An. We then restrict the results of two loop calculations in 0, as given in (A.l), to 
a supersymmetric theory with ny vector superfields and ns chiral superfields. The metric 
and one-form become 

ds 2 = g IJ (g)dg I dg J = 4n v ^(dg) 2 (l + og 2 ) + § dY^ k dY ijk , 

9 (5) 
W I (g)dg I = 2n v ~d$(l + \ag 2 ) + ±(Y^ k dY ljk + dY^ k Y ijk ) . 

The parameter a, which is linear in R and C (where tr(t a tb) = —RS a b and C is similarly 
defined for the adjoint representation, C = N for G = SU(N)), is scheme dependent. To 
the order given in ([5]) the metric is clearly flat and the one form is exact. The leading 
perturbative result is also manifestly positive which should therefore remain valid at least 
in some region of weak coupling. The corresponding 3-loop order result (A. 5) for a, which 
is independent of a because of cancellation between the the a and Wiff 1 terms of (^), 
becomes 

8a = 8a + 2n^C(3C-^ 4 -4tr(t 2 t 2 #-tr(t 2 (yf))^ 2 , (YYfj = Y Um Y jlm , (6) 

with ao the result for free fields, 8ao = l(9ray + ns)- For asymptotically free theories 
ao = cluVi the value at the UV fixed point. Applying (H) to pure SSYM theory without 
chiral matter gives 8a = |nv r (l+4C 2 ^ 4 ) which increases away from the free UV fixed point. 
This shows that a itself does not define a good c-function describing monotonic RG flow 
between fixed points. The three loop corrections in both (A. 5) and (|6|) are proportional to 
the two loop terms in the gauge beta function and are scheme independent. 

In this paper we require expressions for the gauge beta function to three loops and for 
the Yukawa beta function to two loops. At this order they are dependent on the choice of 
renormalisation scheme. For the gauge /3-function we use the NSVZ form, which expresses 
/3 9 to all orders in terms of the anomalous dimension matrix 7^ for the chiral fields and 



has been verified at three and four loops in [12]. This is here written as 



Vfo)= ^^3 , /3(£) = -(/3 -2tF( 7 t 2 ))£ 2 , /3 = 3C- J R, (7) 

where tr = n^tr. The Yukawa beta function is also in general directly expressible in terms 
of 7, as a consequence of the non-renormalisation of the superpotential, 

(3 Y ijk =3Y e(ijl * k) , f l3k = 3^eYiW . (8) 
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To one and two loop order using dimensional reduction [p~3| B 



1 a)i j = pi j = i { YYy j + 2g\t*y j , 

± (9) 
7 (2)< . = _ Y ike Y jkm P m e + 2g 2 (Pt 2 Y J + 2(3 g 4 (t 2 Y J . 

From (A. 9) in the same dimensional reduction scheme 

a = -2C + 5/3 . (10) 

We first show that (|6|) at a perturbative fixed point agrees with to 3-loop order. 
Using Y-(3 Y = 3tr((Fy)7) and the one loop form for 7 in (|9|) it is easy to see that 

8a = 8a - tr( 77 ) + 2n v C(3 g A + ± Y-(3* . (11) 

To obtain a perturbative fixed point it is necessary to choose the gauge group G and 
representation content of other fields so that (3q is effectively small. This is possible |l4j] in 
a large N limit, where as N — > 00, the gauge group is such that C = O(N), ny = 0(iV 2 ), 
and (3q > is tuned to be 0(1). In this case there is a IR fixed point with g 2 = 0(N~ 2 ) 
and the other couplings are given in terms of g# so that in the case of N = 1 SUSY theories 
discussed above Y*,Y* = 0(g*)- I n these circumstances 7 = 0(iV -1 ) and tr( 77 ) = 0(1) 
while the last two terms in (|TT|) are of higher order in 1/N. It is evident therefore that 
( pTTj) is compatible with (@) to leading order. Equivalently in discussing a perturbative 
fixed point it is legitimate to assume that the lowest order beta functions which actually 
appear in ( |TTD are 'completed' by higher order terms, —floQ 2 ~^ P(g)(l + 0(g 2 )), which 



make them vanish. This allows the last two terms in (11) to be dropped and we then 
observe agreement with (0) to 3-loop order. 

We now demonstrate that the formalism can be extended in a fairly simple way to 
obtain results for a to four-loop order, which can then be compared with the exact formula 
(H) . The main idea is to use the integrability conditions for (|3|) to generate the three-loop 
contributions to gu and Wi in the scheme in which the beta functions are presented above. 
We begin with the the g Y y sector and then extend results to the complete metric. The 



1 The compatibility of these results with the NSVZ /3-function with independent 3 loop 



calculations for /3 9 , up to the freedom corresponding to a change of scheme, was shown in [12,13]. 
The form assumed (||) does not itself remove possible scheme dependence since taking SYijk = 
3Y e(ij h £ k) , 5Y ijk = 3h^ £ Y jk)e and using the general result for the variation of the /3-function, 
Sp 1 = ^djSg 1 - SgJdjfi 1 , gives in (§) 5 7 = 0-dh - Sg-dj + [~/,h]. At a fixed point then 
£7* = [7*,/),*] so that the eigenvalues of 7* are invariant. However this does not maintain the 
NSVZ form in (|) which appears to define the renormalisation scheme uniquely. 
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Fig. 1 Two, three and four loop vacuum diagrams for chiral scalar fields. Each 
diagram represents a potential contribution to a, although the results here contain 
no term corresponding to the non planar diagram, while selecting one vertex de- 
termines the possible forms for the one form W and pairs of vertices the metric g. 
Since there is just a single diagram at two and three loops W is exact to this order 
while there are three possible terms in the metric at three loops. 

form of possible contributions are determined by vacuum diagrams, exhibited in Fig. 1 for 
the cubic interactions determined by the superpotential. Hence, for suitable coefficients 
a, P, 7, e (which in general are scheme dependent!) , it may be written as 



dY-g YY -dY = ldY ijk dY ljk + atv((dYdY)(YY)) + /3tr((dff)(fdf)) 

+ 7 tr((YdY)(YdY) + (dYY)(dYY)) + eghv ((dYdY)t 2 ) . 



(12) 



In a similar fashion we may extend the result in (|5]) for Wy-dY + Wy'dY to allow for 3 
loop contributions. However it remains exact to this order and may be disregarded so that 
t YY ~~ ^yy ~~ QyY' -^-t the same effective order we may also take 



dY-t Yg dg + dY-t Yg dg = 2 V gdgtr(d(YY)t 2 ) . 



(13) 



Using (|T^) and (|T3|) with the the results for (3 Y provided by (|8|) and (|9|) then gives an 
equation determining the dependence of a on f, 

dY-d Y 8a = ±tr((dYY)(YY)) + 2gHr((dYY)t 2 ) 

+ \adY M Y jkm {YY) m ,{YYy % + \{j3 + 2 7 - l)Y M Y jkm (dYY) m e (YYy % 
+ \{a + (3 + 2 7 )tr((dyT)(yr)(yT)) + (2a + ^dY^Y^iYY)^ 2 )^ 
+ 4ag 2 Y iki Y jkm (dYY) m e (t 2 y \ + (3a + \e + 2)g 2 tr((dYY)(YY)t 2 ) 
+ 2eg 4 dY M Y jkm (t 2 ) m i(t 2 y i + (e + 4)g 4 tr((dYY)t 2 t 2 ) 

+ 2(l-r 1 )(3 ghr((dYY)t 2 ). (14) 

The first line here follows just from (^) and the one loop result in (|^) . Integrability of flT3| ) 
imposes the conditions 

/3 + 2 7 -l = 2a, e = 4a, (15) 



If, in terms of the previous footnote, we take h = X(YY) then 5a = 8(5 = — 4A, <5 7 = — 2A, 

j 



while <5 7 (2)l i = 4A£ 2 ((YY)t 2 )%- + 8Xg 2 Y M Y jkm (t 2 ) m e 
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and then (I|) then gives 



8a = 8a + \tv((YY)(YY)) + 2gHr((YY)t 2 ) (l + (1 - r])^ 2 ) 

+ 2a Y iU Y jkm P m iPU + ±(3a + l)tv ((YY)(YY)(YY)) (16) 
+ (2a + l)g 2 tr((YY)(YY)t 2 ) + 4(a + l)# 4 tr((yf)t 2 t 2 ) + . . . . 

It remains to consider also the dependence on the gauge coupling. To this end we 
extend the metric in (|5|) to the form 

g gg dg 2 = An v ^(dg) 2 (l + ag 2 + rg 2 ^{(YY)t 2 ) + 0(£ 4 )) , (17) 

as well as requiring 

dgt g ydY + dgt g Y-dY = 2Kgdgtr(d(YY)t 2 ) . (18) 
Then (|), with (0), gives 



(19) 



gd g 8a = 4n v p(g) (l + (2C + a)g 2 + Tg 2 tv((YY)t 2 ) + 0(g 4 )) 
+ 2k g 2 tv((p f Y)t 2 + (Y(3*)t 2 ) . 

It is of course crucial that the Y, Y dependent terms obtained from integrating (|19[) are 
consistent with those already found (|16|). This requires, after using (fLOP, 



«-2a = l, t-2?] = 3. (20) 
We may then calculate the purely g dependent terms in a from QI§|) to be 

8a = -2n v (3 Q g 2 (l + §/3 £ 2 ) + 4gHr(t 2 t 2 ) + f g 6 tr(t 2 t 2 t 2 ) + O((3 g 6 ) . (21) 
Combining (|T6| ) and ( pl|) we can express a to effectively 4 loop order in the form 



8a = 8a - 2n v ^g 2 {\ + §/3 # 2 ) + tr(PP) 

+ 2a Y ike Y jkm P m e P j i + § (3a + 1) tr(PPP) - 4a g 2 tr(PPt 2 ) (22) 
+ 0(/^ 6 ) + 0(/^ 4 )tr(Pt 2 ) . 

For finite theories with zero beta functions a = ao is independent of the couplings and this 
is reflected in ([22] ) since the corrections vanish if (3q = 0, P = which are sufficient for a 
finite theory to two loop order. 
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To make the connection with the results in || we use, to the required order, 



n v P(g) + \ Y-f3 Y = -n v (3 g 2 + tr(PP) - Y M Y jkm P m e P\ 

+ 2ghr(PPt 2 ) + 2l3 ghi{Pt 2 ) , (23) 

tr( 77 ) = tr(PP) - 2Y iki Y jkm P m e P^ + 4g 2 tr(PPt 2 ) , 



and hence (|22|) becomes consistent with 



8d = 8a - tr( 77 ) + | tr( 777 ) 



+ 2n v ~m (1 + O((3 g 2 )) + | Y-f + a tr((F/3 f )P) 



(24) 



At a fixed point the second line, and hence the dependence on the undetermined parameter 
a, vanishes. The result (|24]) is in accord with @ taking 7 at the IR fixed point to have 



eigenvalues 7r . It is a tribute to the consistency of quantum field theory that two such 
different approaches give identical results — the non-perturbative method of 0] effectively 
uses only 1-loop triangle anomalies while the second method strictly requires 4-loop cal- 
culations in curved space which were bypassed here by use of the consistency conditions 
associated with @. 

As a specific illustration we consider an example based on magnetic SQCD 0, as 
described by Kogan et al [p~5[| , which for G = SU (N) contains chiral fields ip l , (pi belonging 
to N, N representations of SU(N), with i = 1, . . . Nf as well as a SU(N) singlet M l j. 
There are then Yukawa interactions preserving SU(Nf) x SU(Nf) symmetry represented 
by the superpotential y0iM i j(p^ and the matrix P defined in (|[|) has the form, in a basis 
for the fields given by (M, (p, (p), 



P 



/Ny 2 l N 2 \ 

' {N f y 2 - ^g 2 )l NNf 
V ( Nf f-N^lf) lNNf ) 



(25) 



For this theory (3q = 3N — Nf and the lowest order formula for a becomes 
8a = ±(9(iV 2 - 1) +2NN f ) 

AT. 9 (26) 

- 2(N 2 - 1)(3N - N f )g 2 + N 2 N f 2 y 4 + 2^j-(NN f y 2 - (N 2 - l)g 2 ) . 



This satisfies, to lowest order since from (|[) ds 2 = 4(A^ 2 — l)dg 2 /g 2 + 4NNfdy 2 , 

g— 8a = 4(iV 2 - 1)0? % y) , —8a = WN f 2 fi{g, y) , (27) 



7 



and has a minimum at the IR fixed point 



NN fV* = N 2 +2N f (iV2 " = 3N ~ N f- ( 28 ) 



It is easy to show that then 



N + 2N f 



S(a uv - a Ia ) = tr(P^) = 3 (3N - N f f . (29) 

Of course these perturbative results are only valid in a large iV limit with Nf tuned to 
ensure SN - N f = 0(1). 

We may also note the restriction to N = 2 SUSY theories which may be defined 
by taking the superpotential to be ^/2grj a x T T a ^ where 77, £ and x are chiral superfields 
transforming according to the adjoint, R and R* representations of the gauge group. In 
this case there is a single coupling g with the perturbative beta function non zero only at 
one loop with (3g = 2(C — R) where tr(T a 2],) = — R5 a },. From (|5]) and ([10]) we have 



g gg dg 2 = Any ^(dg) 2 (l + A(3 g 2 ) , (30) 

and, with the analogous result for W, a = cio to three loops, independent of g. The 
finiteness of N = 2 theories beyond one loop, in particular the vanishing of 7, suggest that 
this may be true to all orders, compatible with there being no perturbative fixed point. 



Appendix A. 

In J7| results were described for the non supersymmetric theories. In order to show 
the connection with above results for SUSY models we describe briefly the necessary tran- 
scription code. For real scalars cj) 1 and n$ Majorana fermions if}, ^ = ifj T C, C = — C T , 
the essential dimensionless couplings for a general renormalisable four dimensional theory 
are g 1 = (g,Ti, Xijkt), with g the gauge coupling, a Yukawa interaction ^T^cf) 1 , where 
Tj is a matrix linear in 75 satisfying CTj = — (CTj) T (fj is also defined by 7 M Fj = 1^7^), 
and a quartic scalar interaction jiXijMftft </> (/> ■ The gauge interactions are specified by 
the gauge group generators tf, a = 1, acting on 0, tp satisfying t% = — tf T , 
CtfC -1 = —tf T (where 7 M t„ = ttlv)- For sucri couplings the metric calculated in |7| can 
be written as 

As 2 = Any ±(dg) 2 (l + A + \ tr(df<dr 4 ) + £ d\ ijM d\ ijU . (A.l) 
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The associated one form to the same order was also 



W I (g)dg I = 2n v Ug{l+\A J^J +^ £ tr(f ^0+^^ ^ Ay M dAy M . (A.2) 

The trace in the terms involving T, and tfP, also involves a sum over spinor indices. Using 
the results for the /3- functions 



^(r^r, + ±(1^ + e£* a )r« + r i ±(t j r j + e^ 2 ) + itr(r i f j )r 

^° = 3 (HC — 2i?^ — \R<p) , -R</> = — tr(t^ 2 ) , = — | tr(t^ 2 ) , 



(A.3) 



# " u 16tt 2 (16tt 2 )< 

where 

/3n = h(llC-2R. h - hRA , R^ = -Ti(t* 2 ), R,,, = -A tFfM , 

(A.4) 



ft = I<7(34<7- lOify -ify) - |tr(^ 2 ^ 2 ) -2tr(^V 2 ) , 

then a can be found, from @ and the results for free fields, to three loops, independent 
of A, to be 

8a = i (l24n y + llr^ + 2n ) + * (r^fta 4 - J tr(^ 2 f,r i ) o 2 ) . (A.5) 

Completing the calculation in to include scalar fields, using dimensional regularisation, 
gives 

^■dim.reg. = 17C — ^Rip — ^ • (A-6) 

The reduction to a general SUSY theory, as described above when = ny+ns, = 
2ns, m ay be obtained by taking 

with ta d the adjoint representation generators and P± = A(l ± 75), and also 



Thus R^p = C + R, R (j} = 2R and (|A~7T|) , (|AT2]) reduce to (EH as well as (|A~5|) leading to 



To apply (|AT6|) we transform to dimensional reduction, as appropriate for supersymmetric 
perturbative calculations, using |16| 167i 2 /g 2 eg = 16iv 2 /g 2 ed + gC which gives 

^4-dim.red. = ^4-dim.reg. — \C = g(39C — 17i?) , O = A dimire d. + |i? = 13C — 5i? . (A. 9) 
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